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Abstract. The Bardeen-Cooper-Schrieffer (BCS) functional has recently re- 
ceived renewed attention as a description of fermionic gases interacting with 
local pairwise interactions. We present here a rigorous analysis of the BCS 
functional for general pair interaction potentials. For both zero and positive 
temperature, we show that the existence of a non-trivial solution of the non- 
linear BCS gap equation is equivalent to the existence of a negative eigenvalue 
of a certain linear operator. From this we conclude the existence of a criti- 
cal temperature below which the BCS pairing wave function does not vanish 
identically. For attractive potentials, we prove that the critical temperature is 
non-zero and exponentially small in the strength of the potential. 



Over the last few years, significant attention has been devoted to ultra-cold gases 
of fermionic atoms, both from the experimental and the theoretical point of view. 
Because of the ability to magnetically tune Feshbach resonances it has become 
possible to vary interatomic interactions via external magnetic fields, and hence 
to control the strength of the interactions among the atoms. This method was 
used in remarkable experiments to study the crossover from the regime of Bose- 
Einstein condensation (BEC) of tightly bound diatomic molecules to the Bardeen- 
Cooper-Schrieffer (BCS) regime of weakly bound Cooper pairs in ultra-cold gases 
of fermionic atoms. For a recent review on this topic we refer to [3 [7]. 

The main purpose of this paper is to give a mathematically precise study of the 
BCS regime. This regime is usually described by the BCS functional, derived by 
Leggett [TD] based on the original work of Bardeen-Cooper-Schrieffer [3] . We note 
that while in the original BCS model simple phonon-induced non-local interaction 
potentials were considered, the interaction potentials are local for atomic Fermi 
gases. In order to allow for a wide range of applications, we will consider here the 
BCS functional for general pair interaction potentials V . 

One purpose of our paper is to obtain necessary and sufficient conditions on V 
for the system to display superfluid behavior. Recall that in the BCS model su- 
perfluidity (or, rather, superconductivity in the original BCS case, where electrons 
in a crystal structure have been considered) is related to the non-vanishing of the 
pairing wave function, describing particle pairs with opposite spin and zero total 
momentum. That is, superfluidity occurs if the paired state is energetically favored 
over the normal state described by the Fermi-Dirac distribution. 

More precisely, for arbitrary temperature T — 1/(3 > a system is in a superfluid 
state if there exists a non-trivial solution A to the BCS gap equation 
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with E(p) — J (p 2 — fj,) 2 + \A(p)\ 2 . Here, /i denotes the chemical potential. Notice 
that the gap equation is highly nonlinear. In TheoremEJbelow we will show that the 
existence of a non-trivial solution to the gap equation is equivalent to the existence 
of a negative eigenvalue of the corresponding linear operator 

Kp^ + V , i^ = ( P ^ M) ____. (1 . 2) 

The operator Kp^ is understood as a multiplication operator in momentum space, 
i.e., p 2 equals —A in configuration space. Observe that in the limit T — > this 
operator reduces to the Schrodinger type operator 

\p*-p\ + V. (1.3) 

This linear characterization of the existence of solutions to the nonlinear BCS equa- 
tion represents a considerable simplification of the analysis. In particular, it proves 
the existence of solutions to for a wide range of interaction potentials V, and 
hence generalizes previous results [U EH E3 EH| valid only for non-local V's under 
suitable assumptions. 

Moreover, the monotonicity of Kp^ in (3 guarantees the existence of a critical 
temperature T Cl with < T c < oo, such that superfluidity occurs whenever T < T c 
and normal behavior for T > T c . In particular, there is a phase transition at T c . 
For positive values of /i, we shall show that T c is strictly positive whenever the 
negative part of V, denoted by V- , is non- vanishing and the positive part V+ is 
small enough in a suitable sense. Moreover, we prove that T c is exponentially small 
in the strength of the interaction potential. More precisely, if the potential is given 
by XV for some fixed V G L 1 nL 3/2 , then T c < e~ 1/A as A -> 0. This generalizes the 
corresponding result in [3, Eq. (3.29)] to a very large class of interaction potentials. 
For negative potentials, also a corresponding lower bound will be shown to hold. In 
a forthcoming work j.9], the precise asymptotics of T c for A — > will be investigated. 

At zero temperature, the existence of a non-trivial solution of the BCS gap 
equation is also related to the non-vanishing of the energy gap between the ground 
state and the first excited state in the effective quadratic BCS Hamiltonian on Fock 
space. For very specific attractive interaction potentials, it has been argued that 
this gap is non- vanishing, see e.g. EH1EI1EH]- ^ e n °t e J however, that the existence 
of such an energy gap is not necessarily implied by the non- vanishing of the solution 
to the BCS equation. In fact, the presence of a gap is not necessary for our analysis. 
We consider general potentials where such an energy gap need not be present, a 
priori. We do prove that the gap is non-vanishing for a certain class of interaction 
potentials, however. 

We note that our results are relevant for actual experiments on cold atomic gases 
as long as the description via the BCS functional is applicable. There is an extensive 
literature on this subject, see e.g. El El EH EH ESI El]- According to Leggett 
[TP] this is the case for weak interactions as long as the range of the potential V 
is much smaller than /i -1 / 2 (in appropriate units), i.e., for weakly coupled dilute 
gases. In particular, "weakly coupled" means that the corresponding Schrodinger 
operator p 2 + V typically does not allow for a negative energy bound state since 
otherwise the system is in the BEC regime of tightly bound bosonic molecules. We 
refer to EH] for a detailed discussion on this question. 

2. Model and Main Results 

The precise definition of the BCS functional considered in this paper is given as 
follows. For the convenience of reader, we describe in the appendix the motivation 
and physical background; the discussion there follows Leggett's derivation in [TP] , 
See also [2] for a detailed study of the BCS approximation to the Hubbard model. 
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We use the standard convention a{p) = (2?r) 3 / 2 L 3 a{x)e lpx dx for the Fourier 
transform. 

DEFINITION 1. Let V denote the set of pairs of functions (7, a), with 7 e 
L^R 3 , (l+p 2 )dp), < 7(p) < 1, and a £ iJ^R 3 , do;), satisfying \a{p)\ 2 < j(p)(l- 
7(p)). Let V S i 3 / 2 (R 3 ,da;) be real-valued and n E R. For T = 1//3 > and 
(7, a) 6 2?, the energy functional Tp is defined as 

^(7, a) = J (p 2 - nh(p)dp + J \a(x)\ 2 V(x)dx-^S( 1 ,a), (2.1) 

where 

5(7, a) =- / [*(p) lns(p) + (l-*(p)) In (l-«(p))] dp, 
with s(p) determined by s(l — s) = 7(1 — 7) — |a| 2 . 

As explained in the Appendix, (|2.ip has the physical interpretation of — i(27r) 3 / 2 
times the pressure of a system of spin 1/2 fermions at temperature T interacting 
via a pair potential given by 2V. 

It is natural to require V e i 3/2 (R 3 ), as we do here, since this guarantees that 
V is relatively form-bounded with respect to p 2 (see, e.g., Sect. 11 in [H]). In fact, 
it would be enough to require the negative part of V to be in I?! 2 + Lf , and allow 
for more general positive parts of V. For simplicity of the presentation, we will not 
pursue this generalization here. 

In the absence of a potential V, it is easy to see that J-p is minimized on T> by 
the choice a = and ^(p) = 70 (p) = [e^ 2- ^ + We refer to this state as the 
normal state of the system. We will be concerned with the question of whether a 
non-vanishing a can lower the energy for given potential V. 

Note that at T — it is natural to formulate the model in terms of a functional 
that depends only on a and not on 7. Namely, for fixed |d(p)| 2 , the optimal choice 
of 7(p) is 

y(0) f 2(1 + Vi-MHp)\ 2 ) for P 2 < H , , 

1{P) ~ \ |(1 - yi^W) for P 2 > » { ' 

in this case. Similarly, also at T > it is enough to minimize J-p only over pairs 
(7, a) satisfying a certain relation independent of V, which is displayed in Eq. (|3.4p 
below. 

According to the usual interpretation of the BCS theory, the normal state is 
unstable if the energy can be lowered by the formation of Cooper pairs, i.e., if Tp 
can be lowered by choosing a ^ 0. The following Theorem, which is the main result 
of this paper, shows that this property is equivalent to the existence of a negative 
eigenvalue of a certain linear operator. 

THEOREM 1. Let V G L 3 / 2 (R 3 ) ; (i <E R, and 00 > T = I > 0. Then the 
following statements are equivalent: 

(i) The normal state (70,0) is unstable under pair formation, i.e., 

inf Fp{j, a) < Tp (70, 0) . 

(7,a)£l5 

(ii) There exists a pair (7, a) € T>, with a ^ 0, such that 

a (p) = -TT^*® ( 2 - 3 ) 
lip) ~ 5 

satisfies the BCS gap equation 

A = -V * (§ tanh ^ ) . with E (P) = Vip 2 -P) 2 + |A(p)| 2 . (2.4) 
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(iii) The linear operator 

/3(p 2 - M ) i 1 

K^ + V , Kfi^tf-ti jv-A-V (2 ' 5) 
has at least one negative eigenvalue. 

Here, * denotes convolution in the BCS gap equation (|2.4|) . In order to avoid 
having to write factors of 2n, we find it convenient to define convolution with a 
factor (2^)~ 3 / 2 in front, i.e., (/ * g)(p) := (2tt)- 3 / 2 / r , f(p - q)g{q)dq. 

Remark 1. The operator (|2.5p is obtained by taking the second derivative of the 
functional Tp with respect to a at the normal state (70,0) (see Eq. ()3.10p below). 
Since the normal state is a critical point of Tp, it is then easy to prove that (iii) 
implies (i). The opposite direction is not immediate, however. It says that if (70, 0) 
is not a minimizer, it can also not be locally stable. 

Remark 2. As mentioned above, a minimizer of J-p satisfies the ^-independent 
relation (2.2) at T = 0, and (|3.4|) at T > 0. It is not necessary to assume that these 
equations hold in part (ii) of Theorem [1] however. 

Note that in the zero temperature case T = the gap equation (|2.4[) takes the 
form 

A = -V*± 

(compare with [3l [HI [lOl [T2] ) and the operator in (|2.5p is given by \p 2 — fi\ + V. 

For any potential V we shall prove the existence of a critical temperature < 
T c = 4 < 00 at which a phase transition from a normal to superfluid state takes 

P c 

place. More precisely, we shall show that the gap equation (|2.4[) has a non-trivial 
solution for all T < T c , whereas it has no non-trivial solution if T > T c . This 
follows from the next Theorem, in combination with Theorem [T] 

THEOREM 2. For any V G L 3/2 (R 3 ) there exists a critical temperature < T c = 
4 < 00 such that 

P c 

inf spec( Kp^ + V ) < if T <T C 
Kp^ + V>0 if T>T C . 

With the aid of the Birman-Schwinger principle, it is possible to characterize the 
value of the critical temperature T c , at least if T c > 0. In order to do this, it is 
convenient to decompose V into its positive and negative parts, i.e., V = V+ — V—, 
with V+ > 0, V- > and V-V+ = 0. We shall show in Proposition H that T c > 
satisfies the equation 

1= V 1 ' 2 - l -^V 1/2 . (2.6) 

Kp*# + V + 

Note that since Kp^ > 2/(3, this leads immediately to the general rough upper 
bound 

T c <U V -\\oc- (2-7) 
For negative potentials V and a positive chemical potential fi, we shall show that 
the critical temperature is always positive, i.e., that \p 2 — fj,\ + V has a negative 
energy bound state. We can even say more, namely we can guarantee superfluidity 
for small temperature whenever the positive part V+ is not too large in a suitable 
sense. 

THEOREM 3. Let V E L 3 / 2 (M 3 ) be not identically zero, and let (j, > 0. 

(i) If V < 0, then the corresponding critical temperature is positive, T c = -i- > 
0. 
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(ii) Let V = V+ — V-, and let (3 > (3 C {—V-), where f3 c {—VJ) is the critical 
inverse temperature for the potential —V- . Then the system at temperature 
1/(3 is still in a superfluid state if \\ V+ 2 K^^V^ 2 \\ is small enough. 

In particular, this theorem states that for any negative potential one can add a 
sufficiently small positive part such that the critical temperature T c is still strictly 
positive. 

The next question concerns the magnitude of T c . The rough upper bound given 
in (|2.7|) turns out to be much too rough, in general. In fact, the following Theorem 
shows that T c is always exponentially small in the strength of V. To be able to 
state the bound precisely, let a :— inf t>0 i(e t + l)/((i + 2)(e* — 1)). Numerically, 
a « 0.654. Let further 

f(t) := At r dpp 2 ( , , 1 , 5— ^ ^ . (2.8) 

JK) 2vr 2 i \\pi-l\ + t p 2 + l + t) V ' 

Note that / is a positive and monotone decreasing function, with f(t) ~ ln(l/i) as 
t — ► 0. We obtain the following upper bound on the critical temperature T c . 

THEOREM 4. Let f~ x be the inverse of the function f in (B.tfjj . Let Li > 0, and 

letV = V+-V-, with V £ L 3 / 2 , V- £ L 1 n L?l 2 , and with ± (f) 4/3 ||Vl|| 3/2 < a « 
0.654. Then the critical temperature satisfies the upper bound 

T < M f -x ^-|(|) 4/3 r-ll3/ 2 \ f2q . 

Note that, in particular, /~ 1 (i) ~ e~* for large t. Hence the critical temperature 
is exponentially small in the potential for V- S L 1 n L 3 / 2 . 

For negative potentials, one can show that the bound given in Theorem is 
indeed optimal. More precisely, if the potential is given by AV for some fixed 
negative V £ L 3 / 2 , we shall show in Proposition [5] the lower bound 

T c > 5m, v (A) , 

where g^y is some positive function depending only on /i and V satisfying g^.v{t) ~ 
e -V f f or small t. Combining both bounds then implies that for negative potentials 
V the critical temperature is exponentially small in the coupling parameter. This 
generalizes well-known calculations in the physics literature, see e.g. [3j Eq. (3.29)], 
to a very large class of interaction potentials. Note that, in particular, this result 
implies that it is not possible to calculate the critical temperature via a perturbative 
expansion in V . 

Finally, we comment on the continuity properties of minimizcrs of J- p. In fact, for 
positive temperature T > the minimizing 7 and a can be shown to be continuous 
functions. (We show this in Prop. [3] below.) This may not be surprising, since 
already the minimizer without potential, jo(p) = l/(e^ p + 1), is continuous. 
But in the case of T = the situation is not so clear, since the normal state is a 
step function with jump at p 2 — fx for fx > 0. 

The continuity of a minimizer 7 at T = is closely related to the existence of 
an energy gap, given by 

S:=inf£(p), (2.10) 

p 

with E defined in (|2.4|) . As explained in the Appendix, E(p) has the interpretation 
of the dispersion relation for quasi-particle excitations. The following shows that 
in case V has a strictly negative Fourier transform, 7 is continuous even at T = 
and the gap S is non- vanishing. 
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Proposition 1. Assume that inf spec(|p 2 — fi\ + V) < 0, and that V is strictly 
negative. Let (7, a) EV be a minimizer of T^. Then S > and 7 is a continuous 
function. 

This result implies, in particular, that 2 > for negative V with strictly negative 
Fourier transform. We note that this property need not necessarily be true for all 
potentials V, however. 

We conclude this section with two remarks. 

Remark 3. Our results can be generalized to non-local potentials given by an 
integral kernel V(x, y). Exactly such a non-local potential was used in the original 
paper of BCS [3j who considered V of the form V(x,y) = — Vo(/)(x)(f>(y), with 
4>{p) = 1 if \p 2 — m| < h and otherwise. This models an effective potential arising 
from electron-phonon interactions in crystals. In this case, the minimizer of can 
be evaluated explicitly (see, e.g., [12]). 

In particular, our analysis can be used to generalize previous results [41 IT31 IT5I [TO] 
on the existence of solutions of the gap equation, as already mentioned in the 
Introduction. None of this previous work allows for an interaction kernel V(k, k') 
of the form V{k — k'), however, since such a kernel does not satisfy appropriate LP 
conditions. Such conditions were necessary for utilization of fixed point arguments 
in the previous investigations of this problem. 

Remark 4. In [H1|T5] the BCS-BEC crossover is studied numerically for certain po- 
tentials V that are negative and have negative Fourier transform. These potentials 
satisfy all the assumptions of our results stated above. 

In the following Sections [3] and dj the proofs of the claims of this section will be 
given. 

3. Properties of the Minimizer of Tp 

In this section we shall give the proof of Theorem [1] and Proposition [1] stated 
above. We start by showing the existence of a minimizer of J- p. 

Proposition 2. There exists a minimizer of Tp in T>. 

Proof. We first show that Tp dominates both the L 1 (R 3 ,(1 + p 2 )dp) norm of 7 
and the ii^R 3 , dx) norm of a. Hence any minimizing sequence will be bounded in 
these norms. We have 



^3(7, a) > Ci + -^j (p 2 - fih( P )dp + J \a(x)\ 2 V(x)dx, 

\ |(p 2 -M)7(p)*-^( 7j a)) =-~ J Hl + e--^ 2 ^)dp. 



where 

d = inf 

(7,q)SX> 

Since V G L?l 2 by assumption, it is relatively bounded with respect to —A (in 
the sense of quadratic forms), and hence Ci = inf spec (p 2 /4 + V) is finite. Using 
\&(p)\ 2 < l(p)i we thus have that 

^P 2 j{p)dp+ [ V (x)\a(x)\ 2 dx > C2 [ r y{p)dp. 



4_ 

Using again that |a(p)| 2 < j{p) < 1, it follows that 

^>(7, ") > ~ A + g\\ a \\H^K 3 ,dx) + gll7lUi(K 3 .(l+P 2 )<*P) > (3-1) 

where 

A = -d-J [p 2 /4 - 3/i/4 - 1/4 + C 2 ] _ dp , 



BCS FUNCTIONAL FOR GENERAL PAIR INTERACTIONS 



7 



with [ • ]_ = min{ • , 0} denoting the negative part. 

To show that a minimizer of exists in T>, we pick a minimizing sequence 
(7„,a n ) £ V, with T^(pj n ,a„) < 0. From (|3.1|) we conclude that ||a n |llri < 8A, 
and hence we can find a subsequence that converges weakly to some a 6 H . Since 
V e L 3/2 (M 3 ), this implies that 

lim / V(x)\a n (x)\ 2 dx= [ V(x)\a(x)\ 2 dx 
ra-too J R3 J R3 

dH Thm. 11.4]. 

It remains to show that the remaining part of the functional, 



^j(7,«) = 



{P 2 ~ V>h(p)dp+ 7} s (p) lns(p) + 4(1 - s (p)) l n (! - s (p)) 



dp, 
(3-2) 

is weakly lower semicontinuous. Note that J-g is jointly convex in (7, a), and that 



its domain I? is a convex set. We already know that a n — 1 a weakly in H 1 ( 
Moreover, since 7„ is uniformly bounded in L 1 (R 3 ) nL°°(R 3 ), we can find a subse- 
quence such that 7„ — 7 weakly in L P (M. 3 ) for some 1 < p < 00. We can then apply 
Mazur's theorem [111 Theorem 2.13] to construct a new sequence as convex combi- 
nations of the old one, which now converges strongly to (7, a) in L P (M 3 ) x L 2 (R 3 ). 
By going to a subsequence, we can also assume that 7„ — > 7 and a n — ► a point- 
wise [HI Theorem 2.7]. Because of convexity of T^, this new sequence is again a 
minimizing sequence. 

Note that the integrand in (13. 2ft is bounded from below independently of 7 and 
a by — j3~ x ln(l + e -/3 ( p _M '). Since this function is integrable, we can apply Fatou's 
Lemma [HJ Lemma 1.7], together with the pointwise convergence, to conclude that 
liminf ^(7„,a„) > ^§(7,5). 

We have thus shown that 

Fa < lim inf Tq (j n , a») • 

n — >oo 

It is easy to see that (7, a) G V, hence it is a minimizer. This proves the claim. □ 

Next, we establish a few auxiliary results that will be needed below in the proof 
of Theorem [TJ 

Lemma 1. Let (7,0) 6 T> be a minimizer of J-p for < j3 < 00. Then, for a.e. 
p e R 3 , 

(V,a)(p) = (p 2 -^^- Y , (3.3) 

(P 2 M) = - 27(P j -1 / (2^( 7 (p)-i)2 + |Q(p)p) , (3.4) 
wftere /(a) := ± In i±2 / or < a < 1. 

Proof. For e > 0, let A e := {pel 3 : |d(p)| 2 + ( 7 (p) - 1/2) 2 < 1/4 - e}. Let 
g e i/ 1 (M 3 ) DL 1 (R 3 ), with Fourier transform supported in A e . Then (7, a+t<?) G P 
for small enough By the Lebesgue dominated convergence theorem, it follows 
that 



d , . 
—Fp('j,a + tg) 



t=o 



= 2 Re y y(a;)g(x)a(a;)da; 



I Re 



ff(p)d(p)/ (2^(7^) - ±) 2 + |<%)| 2 ) dp. (3.5) 
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d . . 
—F /3 {'y + tg,a) 



t=o 



Similarly, for g real, (7 + tg, a) £ T> for small \t\, and we get 

(P 2 - ^)g(p)dp (3.6) 

+ 1 Jg(p)(2j(p) 1)/ f2^/( 7 (p)-i)2 + |d(p)|^ dp. 

Since (7, a) minimize by assumption, the expressions in (|3.5p and (|3.6p vanish. 
It follows that for a.e. p e i e we have 

(V- * a)(p) = -ia(p)/(2 v /( 7 (p)-i) 2 + |a(p)| 2 ) , 

(P 2 - M) = "^(27(P) - l)/(2 v /( 7 (p)-|) 2 + |a(p)| 2 ) . (3.7) 

To conclude that (13. 7| holds a.e. in M 3 , it remains to show that B := M. 3 \\J e> oA e 
has zero measure. By definition, if p 6 B we have 7 (p)(l — 7(f)) = |<S(p)| 2 i i- e - 
s(p)(s(p) — l) = 0. It is then easy to see that if B has non-zero measure, ^3(7, a) 
can be lowered by modifying 7 and aonB (since 1 1— > i In i has a diverging derivative 
at zero). This contradicts the fact the (7, a) is a minimizer, and hence £? has zero 
measure. □ 

Remark 5. A similar strategy as in the proof of Lemma [T] leads to the conclusion 
that Eq. Q3.3p is also valid at T = 0. The relation between 7 and a, given in 
Eq. (pLU) for T > 0, has to be replaced by ((221) at T = 0, however. We omit the 
details. 

Note that Eq. (|3.4[) and the positivity of / immediately imply that j(p) > 1/2 
for p 2 < fi, and 7(73) < 1/2 for p 2 > Of greater importance below will be the 
following monotonicity property, which bounds 7 in the opposite direction. 

Lemma 2. Assume that (7,0;) <E V satisfies Eq. I[3.4\ ), and let 7o(p) = [e^ p ~^ + 
1] . For a.e. p 6 R 3 , we have lip) < 70 (p) */p 2 < M> for p 2 > pt we have 

lip) > 7o(f0- Moreover, these inequalities are strict if a(p) 7^ 0. In particular, 

. ^ - — y-r < ~ ~ — z — r 071 the Support of OL. 

l-2 70 (p) l-2 7 (p) 

Proof. Let p 2 > /i such that cii(p) ^ and assume by way of contradiction that 
7(p) < 7o(p)- We then have 



\ 7o(p) < v /|A(p)| 2 + (i-7b)) 2 . 



2 

Since / is a strictly monotone increasing function, we obtain 



fi\ - 2 70 (P)) < / ^|a(p)| 2 + (i- 7 (p)) 2 
and hence 

(1 - 2to(p))/(1 - 2 7o (p)) < (1 - 2 7 (p))/ ^2 v /|a(p)|2 + (I- 7 (p))2^ . (3.8) 

This, however, contradicts the fact that both (70,0) and (7, a) are solutions of 
(13. 4|) . and hence both sides of (|3.8p are equal to /3(p 2 — 

The case p 2 < /i is treated similarly. □ 

We now have all the necessary prerequisites at our disposal to give the proof of 
Theorem [TJ 
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Proof of Theorem [JJ [i] =>• [ii] : Consider first the case T > 0. According to Lemma[T] 
a minimizing pair (7, a) £ V with a ^ satisfies Eqs. ()3.3|) and (|3.4p . Using the 
definitions of A and E in $n$) and we see that >/(7(p) - 1/2) 2 + |5(p)| 2 = 
1/2 ~J {p) E(p). Plugging this into Eq. (JO} leads to 



27 (p) - 1 /T V - A* 7 /9^(p) 2 70) ~ 1 V P 2 - M 



or 



tMi « = i^lW E(p) . (3.9) 
2 p 2 - — p, 



Therefore, we can rewrite the relation between 5 and A as 5(p) — 2E(p) t ann ■ 
Together with Eq. (|3.3[) this gives the required form of the gap equation. 
The case T = is treated similarly (cf. Remark [5]). 

[ii] [iii]: Given T > and a A satisfying (|2.3p for some (7, a) £ T>, we 
first show that there exists a (7, 5) £ T> yielding the same A but satisfying in 
addition Eq. (|2.2| at T = or Eq. (|3.4p at T > 0, respectively. Namely, we can let 
a(j>) = m(p)a(p) and j(p) = 1/2 + m(p)(j(p) — 1/2) without changing A. If we 
choose 

t ^ ( P—A_\ tanh(/3£;(p)/2) 
m{p) = sgn — - — - , 

v 1 -^); 2 V /( 7 ( P )-i) 2 + id(p)|2 

with -E(p) as in (|2.4[) , it is easy to see that (|3.9p is satisfied. Using in addition the 
gap equation, the same reasoning as in the previous part of the proof leads to the 
conclusion that Eqs. ()3.3p and (|2.2p (respectively Q3.4p ) are satisfied for this (7,5). 

Hence there exists a (7, a) £ V satisfying (|3.3p as well as (|2.2p (at T = 0) or 
(|3.4p (at T > 0). Using (|3.3p and the definition of Kp,/M, we have that 

1 1 

l-2 l0 (p) ~ l-2 7 (p), 
It follows from Lemma [2] that this expression is negative whenever a is not identi- 
cally zero. 

[iii] =>■ [i]: Consider first the case T > 0. Showing that the existence of a negative 
eigenvalue of the operator (|2.5p implies that the minimizing a is not identically zero 
is equivalent to proving that a = implies that inf spec (Kp tfi + V) > 0. First, we 
note that, for any 8 £ C£°, \\(t8) 2 + (70 - ;j) 2 ||oo < \ for all small \t\. With the aid 
of the Lebesgue dominated convergence theorem it is easy to see that ^2^/3(70, to) 
exists for small t, and is given by 

d 2 

■^^0(70, ts) 

= 2 1 V{x)\8{x)\ 2 dx 



(a, (K^ + V)a) - < a, ( - ) (j? - y)a 



2 



pj i\ts(p)\ 2 +( 10 (p)-mho(p)(i-jo(p))-\t6(p)\ 2 } p - 

Assuming that the minimizing a is identically zero, it follows that the corresponding 

d 2 

d 2 



7 = 70 and hence 4r ^73(70, t6)\t=o > 0. But 



dt 2^{lo,tS) 



t=0 



= 2 / V(x)\S(x)\ 2 dx + 2 / \8(p)\ 2 K^(p)dp, (3.10) 



and hence (8, {Kp tll + V)8) > for all 8 £ Cq°. This proves the statement. 
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In the case T — 0, it is sufficient to minimize Tp over (7, a) S T> with 7 of the 
form (|2.2p . as remarked earlier. For simplicity, we abuse the notation slightly and 
denote the resulting functional by J-^a). Subtracting its value for a = 0, it is 
given by 

^oo(a)-^oo(0) - \ J \p 2 -fi\(l-^l-4\a(p)\2)dp + J V(x)\a(x)\ 2 dx. (3.11) 

Assume that |p 2 — ^i|+U has a negative eigenvalue. We can then find an a e (7q°(K 3 ) 
such that (a, (\p 2 — [i\ + V)a) < 0. For < e <C 1 small enough, we then have 

^oo(ea) - ^oo(O) = e 2 (a, (\p 2 - f i\ + V)a) + 0(e 4 ) < . (3.12) 

This implies that a minimizer a does not vanish identically. □ 

In the following, we investigate the continuity properties of a minimizing pair 
(7, a) eP. 

Proposition 3. // (7, a) S T> is a minimizer of !Fp for < (3 < 00, then both 7 
and a are continuous functions. 

Proof. Since V S L 3/2 and a e H 1 . we know that V £ L 3 and a £ L 6 / 5+e for all 
e > 0. Hence V *a is continuous fTT] Thm. 2.20]. From Eq. (|3.3|) and the definition 
of A and E given in (|2.3p and (|2.4|) . we observe that y * a — — 5A and hence both 
A and E are continuous functions. As in Eq. (|3.9j) in the proof of Thm. [2 Eq. (|3.4|) 
can be rewritten in the form 

v 2 - a e^P) + 1 

>' " bL -E(p). (3.13) 



2j(p) - 1 el 3 ^) - 1 ' 

The right side of this equation is continuous and does not vanish, not even when 
E{jp) = (which could happen, in principle, for p 2 = /j,). Hence 7 is continuous. 
Using again Eq. (|3.3p , this also implies continuity of a. □ 

The situation is different at zero temperature, however. In the limit f3 — > 00, 
the right side of Eq. (|3.13p vanishes in case E(p) = 0. Hence 7 could possibly 
be discontinuous if p 2 = \i. In fact, this is what happens in the case of vanishing 
potential. On the other hand, a non- vanishing gap S = inf p E(p) > implies 
continuity of 7. At T = 0, the gap equation (|2.4p takes the form 

E(p)a(p) = -(y*a)(p), (3.14) 

\p 2 - mI 

with E(p) = — . Hence the gap 5 is guaranteed to be non-zero if the 

y/l-4\a{p)\ 2 

right side of (|3.14p vanishes nowhere on the sphere p 2 = [i. This can be easily 
shown in case V is strictly negative, which is the content of our Proposition [1] 

Proof of Proposition^ Since V < 0, 

&(V*a) > j \&\(V* \a\). (3.15) 



Moreover, since V is assumed to be strictly negative, there can be equality in (|3.15p 
only if d(p) = e IK |a(p)| for some constant keI. This implies that this is the case 
for a minimizer of J-qc, and hence V * a vanishes nowhere. The statement then 
follows from Eq. (|3~T4|) . □ 
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4. The Critical Temperature 

In this section we investigate the relation between properties of the interaction 
potential V and the critical temperature T c at which the phase transition takes 
place. In particular, we shall prove Theorems [2] 2) 

Proof of Theorem^ Let j3 c = Mp{f3 > : inf spec(Kp^ + V) < 0}. It is easy to 
see that f3 c > 0, since V is relatively form-bounded with respect to \p 2 — /i and 
> c (|p 2 — ^| + 1//3) for some c > (cf. Lemma [3J below). 
From the definition of Kp tll in (|2.5p one easily sees that 

K Pulx <K^ if ft > ft. (4.1) 

It follows that inf spec(Kp^ + V) < for all (3 > ft. Moreover, since Kp^ depends 
continuously on (3 (in a norm resolvent sense), we have that Kp c ^ + V > 0. Setting 
ft — oo in the case inf spec(Kp^ + V) > for all (3 and letting T c = 1/(3 C finishes 
the proof. □ 

Next, we prove the first assertion of Theorem [31 namely that for any negative 
potential the critical temperature is strictly positive. 

Proof of Theorem\3(i) . In light of Theorem [2] it suffices to prove that inf spec (\p 2 — 
/i| + V) < 0. According to the Birman-Schwinger principle, this is equivalent to the 
fact that for some e > 0, the Birman-Schwinger operator 

has an eigenvalue 1. Notice that lim^oo ||-B e || = since V G L 3 / 2 is relatively 
bounded with respect to p 2 (in the sense of quadratic forms). Because of continuity 
in e, it thus suffices to find an e > such that ||-B e || > L 

Hence we have to show that there exists a <f> € L 2 (M. 3 ) with \\<fi\\2 = 1 such 
that (<p, B e (j)) > 1 for an appropriate e. For this purpose, we choose a normalized 

function <f) in such a way that c^yi 1 / 2 does not vanish in a neighborhood of some 
point p with p 2 = /i. We can, for instance, choose 4> to be rapidly decaying, in 

which case l^l^l 1 / 2 ! is continuous and hence strictly positive on some small ball. 
This ball can then be made to overlap with the sphere p 2 = \i by simply multiplying 

<j) by a factor e lqx for appropriate q, which translates the function ^jV"! 1 / 2 by q. In 
this way, 

(4>, B e ct>) = f \<f\V^ 2 ( P )\ 2 1 dp -> oo as e^0. 
J \p l - fi\ + e 

Thus there exists e > such that (</), B e <fi) > 1. □ 

Before proceeding with the proof of the second part of Theorem [31 we introduce 
an alternative characterization for the critical temperature T c in the case T c > 0, 
which we already referred to in Eq. (|2.6p . 

Proposition 4. For all V <E L 3 / 2 (M 3 ), such that V = V + -V-, with V+ > 0, 
V- > and V-V+ — 0, the critical temperature T c — 1/ (3 C is characterized by 



vl /2 - vl /2 



= 1 . 



Kp^ + V+ 
whenever T c =/= 0. 

Recall that Kp^ > 2/(3, hence Kp^ + V+ has a bounded inverse. 
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Proof. Let the operator B@ be defined as 

B@ = V 1/2 - V 1/2 

This operator is compact (in fact, Hilbert-Schmidt) and depends continuously on 
(3 and e. Moreover, \\B^\\ is strictly decreasing in e and 1/(3. 

Since T c > by assumption, we know from Theorem [2] that Kp^ + V has a 
negative eigenvalue for T < T c . Hence the Birman-Schwinger principle implies that 
\\B% || > 1 for 1//3 and e small enough. On the other hand, 1 1 1 1 -> as [3 -> 0. 
Thus there exists a unique /3 such that ||£?g || = 1. 

For (3 > /?, ||-Bf || > 1 for e small enough, and hence Kp tlx + V has a negative 
energy bound state. For (3 < (3, however, \\Bg | < 1 for all e > 0, and hence 
Kp tll + V > 0. This shows that [3 = /3 C . □ 

From Prop. 0] it is then easy to conclude that the critical temperature stays 
strictly positive under small positive perturbations to negative potentials. 

Proof of Theorem\3(ii) . According to Prop. [H it suffices to prove that for (3 > 
(3 C {-V-) and for A > small enough, \\vl /2 (K f3 ^ + XV+^Vl^W > 1 - In °rder to 
simplify the notation let B\ = vl /2 (K f3 ^ + AV+)" 1 K L/2 . It follows that 

u, = r J j ( l - | 1 y 1 / 2 

■: lh I I + A||/v'T:/V + ^ /2 | N 

Since V+ £ I? I 2 by assumption, it is relatively form-bounded with respect to Kp^, 
and hence \\K^J 2 V+K~^ /2 \\ = ||F+ /2 i^7+ /2 || is finite. Therefore ||B A || > 1 if 
1 1 Bo 1 1 > 1 and A is small enough. □ 

Next, we shall derive upper and lower bounds on the critical temperature. These 
bounds are similar to the ones found in the literature, but generalize these results 
from simple rank 1 potentials as in [3] to very general pair interactions. 

We start with the upper bound given in Theorem 21 To this aim we first prove 
two auxiliary lemmas. 

Lemma 3. Let K ^(p) = (p 2 - fx)(e^ 2 -^ + l)(e^ p2 ~^ - l)" 1 . Then there is a 
constant a > such that 

\P 2 ~ Ml + f \ < K , p (j>) < [\p 2 - n\ + . (4.2) 

Proof. Let g(t) = t(e t + l)/((i + 2)(e'-l)) for t > 0. It is easy to see that g(t) < 1. 
Moreover, g(0) = 1 and lim t ^ 00 g(f) = 1, hence a :— inf t> o<?(£) > 0. In fact, 
a w 0.654 numerically. Since Kp^(p) = (\p 2 — fi\ + 2/(3)g((3\p 2 — this implies 
the statement. □ 

Lemma 4. Assume that V G L 1 n L?l 2 and that fj, > 0. Then 

W\ 1/2 , 5 V W\ 1/2 ' 
\P - H\ + e 

where f is defined in (2.8\) . 
Proof. Define k e ^{p) > by 

' ' KM- (4-4) 



<^\\V\\ 1 f{e/ l u) + \ f^j \\V\\ 3/2 , (4.3) 



\p 2 — [i\ + e p ; 



2 i „ i P ^ ne >» 
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Then 
l^l 1 / 2 



\p 2 



-\V\ 



1/2 



< 



\v\ 



1/2. 



1 



P Z + /i 



-\v\ 



1/2 



\v\^K,{ P )\v\^ 



The last term is bounded by 



\V\^k^{p)\V\^ 



(4.5) 



sup 



< (2^)- 3 / 2 ||^|| 00 ||y|| 1 < (2tt)- 

(27r) 3 ^ 1 / 2 /(e/^i), with / given in 



1 1 oo || » ||1 _i \."" J ||"»e,/i||l| 

Moreover, ||fce,nl|i 

For the first term on the right side of (|4.5|) . the use of the Hardy-Littlewood- 
Sobolev inequality and Holder's inequality yields 



\V\ 



1/2. 



1 



p z + fx + e 



\V\ 



1/2 



< sup 

Wl|2=l 



^\v\^,\m 1/2 )<U- 



4/3 



\\V\\3/2 



□ 



Equipped with these two lemmas we are now able to prove the upper bound on 
the critical temperature stated in Theorem UJ 

Proof of Theorem^ Assume that T c > 0, for otherwise there is nothing to prove. 
According to Proposition [4J the critical temperature T c — 1/(3 C is determined by 
the equation 



V 



1/2 



1 



-V_ 



1/2 



= 1 



Using Lemmas [3] and H] as well as V+ > 0, we thus have 



1 < 



V 



1/2 



\p 2 - M | + 2//3, 



■V 



1/2 



1 



4/3 



a 3a \ir 

Since / is monotone decreasing, this implies the statement 



\\V-h/2- 



□ 



Finally, we turn to a lower bound for the critical temperature. We assume that 
the interaction potential is given by XV, with V < not identically zero, and A > 
a coupling parameter. 

Proposition 5. Let V £ L 3 / 2 , with V < not identically zero, and let fi > 0. Then 
there exists a monotone increasing strictly positive function g^y, with g^y(t) ~ 
e -1 /* as t — > 0, such that the critical temperature for the potential XV satisfies the 
lower bound 

T c > ff„,v(A) (4.6) 

for A > 0. 

Proof. We have already shown in Thm. [3] that T c > for positive A. By Prop. HI 
the fact that V is non-positive and (j4.2|) we have that 

1 



m V2 1 | y |l/2 



> 



1A 

For any function ip £ L 2 with HV'lh 

1/A > / iAvf/ 2 (p) 



\v\ 1/2 - t 



\p 2 - m| + 2//3 c 
1 we therefore obtain 

■ dp . 



\V\ 



1/2 



|p 2 - M | + 2//? c 

The latter integral in monotone increasing in fj c and diverges logarithmically as 
(3 C — > oo for appropriate tp such that V'I^I 1//2 (p) does not vanish identically if 
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p 2 = fi. (Sec the proof of Thm. for an argument concerning the existence of 
such a ip.) This implies the statement. □ 
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Appendix A. Motivation and Physical Background 

In this Appendix, we shall briefly explain the derivation of the BCS functional 
defined in Definition [TJ Our presentation follows the one of Leggett in [TUJ. Similar 
derivations can be found in [HI H2] . The BCS functional is obtained from several 
approximations and assumptions on the underlying many-body problem. We as- 
sume the reader to be familiar with the usual Fock space formalism in quantum 
statistical mechanics. 

Consider a system of spin i fermions confined to a cubic box C K 3 , with 
periodic boundary conditions, interacting via a two-body potential V. Assume, for 
simplicity, that V £ L 1 (M 3 ) so that it has a bounded Fourier transform, denoted 
by V(p) = (2tt)- 3 / 2 / V{x)e~ lpx dx. The many-body Hamiltonian on Fock space is 
then given as 

Here, a\, and a k ,a denote the usual creation and annihilation operators for a parti- 
cle of momentum k and spin a, satisfying the canonical anticommutation relations. 
The momentum sums run over k 6 (2-7r/L)Z 3 , where L is the side length of the 
cube fi, whereas the spin sums are over a € {til}- Units are chosen such that 
h = 2m = 1, m denoting the particle mass. 

Following J3[ [TU] we restrict H to BCS type states. These are quasi-free states 
that do not have a fixed number of particles. For details on the formalism, we refer 
the reader to [2]. Quasi- free states are linear maps p from bounded operators on 
the Fock space T to the complex numbers, satisfying p(A^ A) > 0, p(I) = 1, and 

/o(eie 2 e 3 e 4 ) = p(eie 2 ) / o(e 3 e 4 ) - p(eie 3 ) p(e 2 e 4 ) + p(e 1 e i ) p(e 2 e 3 ) , 

whenever is either a,k,a or a\. . 

Assuming both translation invariance and SU(2) invariance for rotations of the 
spins, the state p is completely determined by the quantities 

70) = p(4,t«m) = P( a Li a k,l) . f( k ) = p{a-k^a kA ) = -p(a^ ktl a k ,t) ■ 

All other expectation values of quadratic expressions are zero under these symmetry 
assumptions. Note that < ^(p) < 1 and \ip(p)\ 2 < ~f(p)(l ~ j(— p)) by Schwarz's 
inequality and the canonical anticommutation relations. Since also ip(p) = tp(—p), 
this implies that \<p(p)\ 2 < 7(p)(l — t(p))- 
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The energy expectation value can be conveniently expressed in terms of these 
quantities: 

p(h) = 2 ]T fc 2 7 (fc) + ^ - pMp) 

k ' ' k,p 

+ n^~"^ (0) (^ 7(?,) ) -^^E^m? -%(*>)• ( A - 2 ) 

The first term in the second line is just a constant times the square of the number 
of particles. The last term is the exchange terms. Both these terms will be dropped 
in the sequel. This truncations can be justified on physical grounds for interaction 
potentials V that are of short range such that s-wave scattering is dominant |10j . 

The von-Neumann entropy of the quasi-free state p under consideration here, 
S(p) = — Trjr/olnp, can be conveniently expressed in terms of the 2x2 matrix 

Note that < T{p) < 1 as an operator on C 2 . The entropy is then given by 
(compare with Eq. (2.C.10)]) 

S(p) = -2j2^c* [T(p)lnT(p)} . 
p 

The factor 2 results from the 2 spin components. 

The thermodynamic pressure, V, of a state at a fixed temperature T = A > 
and chemical potential p is defined as 

P(P) = P| [TS(p)- P (M-pN)} . 

Here, N = k a k a ak,a denotes the number operator. Dropping the last two terms 
in (IA.2j) and taking a formal infinite volume limit, one thus obtains 

P( 7) a) := lim V(p) = --^ / (TTr c *T(p) ml» + (p 2 - M ) 7 (p)) dp 
~ 7T"\3 / \ a (x)\ 2 V{x)dx. 

(27Tj J J K 3 

Here, a denotes the inverse Fourier transform of <p, a(x) = (2tt)~ 3 / 2 J ip(p)e tpx dp. 
If we set J- (3 = —^(2tt) 3 V and replace V by 2V to slightly simplify the notation, 
we are thus led to the definition of the BCS model given in Definition [T] 

Finally, we will briefly explain in which sense E(p) is an effective dispersion 
relation and how it is related to an energy gap between the ground state and the 
first excited state of the system. 

With ip = a the Fourier transform of a maximizer of the pressure functional at 
T = 0, we introduce the following quadratic Hamiltonian as an approximation to 
the full Hamiltonian H: 



E fc2 a L a fc,°- + E ^( fc ) (^ fc ~~ p^ a -p^ a pd + <p(p) a l- 



(7, ft: fc,p 

(A.4) 

The Hamiltonian H qu — pN can be diagonalized in a standard way (see, e.g., [T^]) 
via a Bogoliubov transformation. Using the fact that cp satisfies the gap equation 
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E(k)ip(k) = -£(27r) 3 / 2 £ fe , V(k - k')tp(k'), where E{k) = ^(k 2 - fi) 2 + \A(k)\ 2 
and A is determined from ip via ()2.2p - (|2.4p . this leads to 

H qu — [M — inf spec (H qu - fiN) = J2 E ( k ) { C W + 4 d fe) • ( A - 5 ) 

k 

The operators Ck, ct, dk and d) k satisfy the canonical anticommutation relations, 
and both Ck and dk annihilate the BCS ground state. Equation (|A.5|) thus explains 
why E(k) is interpreted as the quasi-particle dispersion relation and 3 = inf^ E(k) 
as the energy gap of the system. 
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